Mean-field glass transition in a model liquid 
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We investigate the liquid-glass phase transition in a system of point-like particles interacting via 
a finite-range attractive potential in D-dimensional space. The phase transition is driven by an 
'entropy crisis' where the available phase space volume collapses dramatically at the transition. 
We describe the general strategy underlying the first-principles replica calculation for this type 
of transition; its application to our model system then allows for an analytic description of the 
liquid-glass phase transition within a mean-field approximation, provided the parameters are chosen 
suitably. We find a transition exhibiting all the features associated with an 'entropy crisis', including 
the characteristic finite jump of the order parameter at the transition while the free energy and its 
first derivative remain continuous. 
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I. INTRODUCTION 

Notable progress in understanding fundamental as- 
pects of structural glasses and their freezing transition 
has been mad* 1 ; significant advancements originate from 
the dynamic formulation of the glass transition using 
mode coupling theory^, while the statistical mechanics 
approach draws extensively from the analogy between the 
well studied spin-glass^ and the glassy solidi While the 
spin glass is characterized by random frozen orientations 
of spins due to the presence of quenched disorder, the 
structural glass is characterized by random frozen space 
positions of the particles but does not rely on the pres- 
ence of quenched disorder. A central idea put forward in 
this context is the concept of the 'entropy crisis— driving 
the transition into the glassy state through a collapse 
of the system's phase space. Transitions of this type 
have shown up in various disordered systems&i and in 
the 'discontinuous' spin-glasses containing no quenched 
disorder 8 (see also Rcf. 9 for similar systems with other 
then spin degrees of freedom). In this paper, we ap- 
ply a heuristic frameworkiS, based on an 'entropy crisis' 
scenario to describe the liquid-glass transition and the 
low-temperature thermodynamics of the glassy state in a 
system of interacting particles in D-dimensional space. 

Progress in our understanding of the liquid-glass phase 
transition and the physics of the low-temperature glass 
state is made along two avenues: i) experimental and 
numerical studies provide new details on specific materi- 
als and on model systems but have little impact on our 
general understanding of the glass phenomenon. On the 
other hand, ii) conceptual studies push our general un- 
derstanding but unfortunately provide little predictive 
power when it comes to the description of realistic sys- 
tems. In choosing a suitable model system, we then have 
to compromise between a realistic description of the glass 
former and one allowing us to make analytical progress, 
e.g., within a mean- field approach. The first model de- 
scribing a liquid-glass phase transition and allowing for 
a mean-field type analysis was proposed by Kirkpatrick 
and Thirumalaiii. Formulated within a density func- 



tional theory, it provided a consistent static and dynamic 
description of the structural glass transition. New in- 
sights into the nature of the glass transition based on 
studies of coupled replicated glassy systemsiS*^ led to 
the formulation of a first principles computational scheme 
providing a description of the equilibrium thermodynam- 
ics of glasses 10 . This scheme has been successfully ap- 
plied, combining analytical and numerical techniques, to 
the soft-spheres model in three dimensions^ - and to the 
Lennard-Jones binary mixture^*. Recently a model with 
point-like particles interacting via a spatially oscillat- 
ing infinite-range potential has been analyzed within this 
framework 15 . Despite its unrealistic structure, this model 
turned out to provide a good testbed for the replica based 
approach: the physics of the liquid-glass phase transition 
revealed in this model turned out fully consistent with 
the heuristic ideas of the 'entropy crisis' scenario. In the 
present paper, we consider a much more realistic model of 
a structural glass: adopting again the scheme suggested 
in Ref. llCt we study the low-density limit of a system of 
particles in -D-dimensional space interacting via a finite- 
range attractive potential of depth Uq, confined to a shell 
of radius R and thickness 2ro, see Fig. ^ It turns out, 
that an appropriate choice of the interaction parameters 
(tq <C R) allows us to adopt a mean-field approximation 
and proceed with an analytic calculation of the free en- 
ergy and a properly defined order parameter to the very 
end. We then arrive at a complete analytic description of 
the liquid-glass phase transition occurring in this model 
glass former and demonstrate that it is again consistent 
with the 'entropy crisis' scenario. 

In spite of this success, we have to admit ignoring vari- 
ous important issues related to the structural glass tran- 
sition: First, we explicitly avoid the discussion of the 
relevance of our equilibrium statistical mechanical ap- 
proach for the apparently non-equilibrium type liquid- 
glass phase transition. This question is both general 
and deep and we are unable to present an answer at the 
present stage. Rather, we hope that, like in the case 
of spin-glasses (another example of matter residing in 
a non-equilibrium state), at least some of the physical 
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FIG. 1: Interaction potential between particles used in our 
model glass former. The attractive interaction (normalized 
to unity) is limited to a shell of width 2ro at a radius R 3> 
ro ^> 1, where we have set the particle radius to unity. The 
repulsive core prevents the collapse of the system into two 
clusters. 



phenomena and observables are well defined and make 
proper physical sense even when computed in terms of 
an equilibrium approach. Second, we also prefer not to 
start a lengthy (and useless) discussion of the question, 
what to do with the crystal state. Such a crystalline 
state is certainly present in our model and, since its en- 
ergy is evidently lower than the typical energy of glassy 
configurations, it is the true ground state of the system. 
In fact, this is the typical situation for most models of 
glasses and the usual algorithm of treating the presence 
of the crystal state (which we also follow in this paper) 
is simple: it has to be ignored. On a qualitative level, 
the reason for such a pragmatic approach is very sim- 
ple: it is well known that frozen glassy states do ex- 
ist at low temperatures regardless of the presence of the 
crystal ground state. Moreover, in many circumstances 
such states turn out to be quite stable during reasonable 
observation times; we then can safely assume that the 
crystal state is located far away from the relevant glassy 
states in configurational space and furthermore that the 
crystal and glassy states are separated by large energy 
barriers. This is a quite standard situation in statisti- 
cal mechanics: except for some rare cases admitting an 
exact solution, one usually studies only a limited ad hoc 
chosen part of the phase space depending on the object 
under study. In the case of the structural glass, one then 
expects two scenaria: either the basic assumption (that 
the existence of the crystal can be ignored) turns out to 
be reasonable and the crystal configurations never show 
up in the calculations. Or, this assumption is wrong and 
then one inevitably faces some kind of instabilities and 
divergences. As for the model and method considered 
in the present paper, the calculations demonstrate that 
the crystal state indeed does not interfere with the glassy 
configurations. 

The paper is structured as follows: In section II be- 
low, we first describe the heuristic framework underlying 
our replica analysis along the lines discussed in Ref. [l(J 
and introduce our specific model. Its free energy is cal- 
culated within the replica mean-field approach in section 



III. The results and conclusions are given in section IV: 
We find that the system freezes at the glass tempera- 
ture T c = Uq/ ln(i?/ro) into an amorphous solid with a 
number of nearest neighbors slightly larger than D. We 
determine the order parameter and its jump at T c and 
present analytical expressions for the free energy and en- 
tropy of the solid and liquid phases as well as for the 
configurational entropy or complexity. 

II. STRUCTURAL GLASS TRANSITION 
A. Symmetry breaking in random systems 

Consider a system of N particles in £>-dimensional 
space described by the Hamiltonian 

1 N 

#M= 2 E ^-^l)' (!) 

where Xj denotes the position of the i-th particle and 
t/(|x|) is the interparticle potential. We assume that at 
low temperatures the system is frozen in a disordered 
(glassy) state which is characterised by random spatial 
positions of the particles; this can be achieved through a 
rapid quench avoiding the crystallization by entropic rea- 
sons. The glass state is characterized by broken trans- 
lational and rotational symmetries, but unlike the or- 
dered crystal configurations (characterized by their spe- 
cific spatial and rotational symmetries), it is impossible 
to identify the residual symmetries left in the glassy state. 
This naturally resembles the spin-glass problem, where 
the spins are frozen in a random state which cannot be 
characterized by any apparent global symmetry break- 
ing. However, unlike spin-glasses, here we do not have 
quenched disorder installed in the initial Hamiltonian. 
Nevertheless, the ideas borrowed from the spin-glass the- 
ory and in particular the use of the replica technique 
turn out to be quite fruitful also for the description of 
structural glasses^. 

In order to demonstrate the effect of spontaneous sym- 
metry breaking, e.g., in ordered magnetic systems, one 
introduces a conjugate field coupled to the order param- 
eter which is set to zero at the end (after taking the 
thermodynamic limit). In spin-glasses, the same strat- 
egy can be applied by introducing several weakly coupled 
copies (replicas) of the original system. Similarly, in or- 
der to demonstrate the freezing of a system of interacting 
particles (described by (JTJ) into a random glass state, we 
introduce two identical copies (with particles at positions 
xj and yi, respectively) of the same system described the 
Hamiltonian 

N N 

H 2 = -£[/(|x l -x,|) + -5>(|y. i -y,|) 

N 

+e^W(x i -y i ). (2) 

i 
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The last term in describes a weak attractive potential 
W between the particles at Xj and y ; of the two systems 
and plays the role of the symmetry breaking conjugate 
field in the ordered system. As usual, the control pa- 
rameter e is set to zero after taking the thermodynamic 
limit and the system can end up in one of two phases: 
i) the particles of the two systems are independent (un- 
correlated), indicating that they do not memorize their 
spatial positions and hence the original system is in the 
high-temperature liquid phase, or ii) the positions of the 
particles remain correlated, indicating that they arc lo- 
calized in space and we conclude that the original system 
is in the low-temperature solid state. As in spin glasses, 
in order to obtain more detailed information about the 
phase transition, it is convenient to introduce to replicas 
of the original system. Also, in the actual calculation 
there is no need to introduce a supplementary attractive 
potential between the replicas: following standard prac- 
tice, it is sufficient to allow for the possibility of symmetry 
breaking in order to prove its existence afterwards. 



B. Entropy crisis scenario 

The entropy crisis scenario^ for the glass transition 
builds on the idea of a phase space collapse upon entering 
the frozen phase; in its pure form it shows up in the ran- 
dom energy model of spin glasses which has been solved 
exactly 6 . We first briefly summarize the main features 
of this heuristic framework as applied to the problem of 
structural glasses, following the original work of Mezard 
and ParisiiS. 

We assume that the partition function Z can be rep- 
resented in the form 



Z = ]Texp(-A7 Q /T), 



(3) 



where Nf a denote the energies of the thermodynami- 
cally relevant local mimima in configurational space. The 
number of local minima with energy / is assumed 

to be exponentially large, fl(f) = exp[NS(f,T)], where 
S(f, T) denotes the configurational entropy density or 
complexity. Finally, the function S(f, T) at fixed T 
shall have the qualitative shape shown in Fig. [21 with 
S(f,T) = at / < f m i n (T), and a concave increase for 

/>/min(T). 

The partition function is written in the form 



Z*t [ dfexp{-N[f-TS(f,T)]/T} 

Jf>fmm 



(4) 



which can be evaluated within a saddle-point approxi- 
mation in the thermodynamic limit (N — ► oo). The free 
energy density F(T) = — T(lnZ) /N is given by 



F(T) = MT)-TS[MT),T], T>T C 



(5) 



with /*(T) defined via d f S\ MT) = 1/T < (3 C and /3 C the 
maximal slope of the function S(f), ci. Fig. [21 At low 




FIG. 2: Qualitative shape of the configurational entropy den- 
sity assumed in the phenomenological description of a glass 
transition driven by an entropy crisis. The slope f3 corre- 
sponds to the inverse temperature and assumes its maximal 
value f} c — 1/T C at the glass transition temperature T c . 



temperatures T < T c the integral is determined by the 
lowest energy state alone, F(T) = / m in(T). 

A microscopic formalism capturing the above phe- 
nomenology can be set up with the help of replicas^: 
Consider to identical (non-coupled) replicas of the same 
system where the particle positions remain correlated 
among the different replicas; this Ansatz describes a 
molecular liquid where each molecule consists of to par- 
ticles originating from different replicas. The partition 
function of the replicated system takes the form 

Z m tt I dfexp{-Nm[f-(T/m)S(f,T)]/T}. (6) 
-'/>/„„„ 

Note that the phase space volume f2(/) remains that of 
the non-replicated liquid since the molecular structure es- 
sentially preserves the configurational degrees of freedom, 
hence the 'entropic temperature' is reduced by a factor 
to. The free energy density F(m,T) = —(T/mN)lnZ m 
of the molecular liquid reads 

r . 7 , ■ f*(m,T) - ^-S[f*(m,T),T], T > mT c 
b (to, T) = < to 



/min(T), 



T < mTr 



where /»(m, T) is defined via the saddle-point equation 



d/S(f,T)\ 



J,(m,T) 



m/T < C . 



(8) 



In order to study the thermodynamic properties of our 
original system (with to = 1) in the low temperature 
glassy phase at T < T c , we continue analytically the 
expression for the free energy density F(m, T) from inte- 
ger values to to arbitrary continuous values and analyze 
its behavior for to < 1. Starting from small to with 
< to < T/T c = to* (T) , the m-replica system resides 
in the molecular liquid phase with a free energy density 
given by © . As to approaches to* (T) from below, the 
free energy density becomes pinned to the lowest value 
/min(T) and the system freezes into the glassy phase. A 
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further increase of m beyond m* (T) results in a constant 
free energy density / m i n (T) and this remains to be the 
case as to — > 1. As a result, the free energy density (and 
hence the entire thermodynamics) of the original glass 
phase can be computed from the m-fold replicated sys- 
tem residing in the molecular liquid phase and taking the 
limit F m=1 {T) = F[m*(T),T]. The critical temperature 
T c is reached when m*(T c ) = 1. At temperatures above 
T c the system resides in a liquid phase: replicas are in- 
dependent, hence m = 1, and the free energy is given by 
F m=1 {T) = F[1,T}. 

The crucial step then is the determination of to*(T). 
Assume we have managed to compute the free energy 
density F(m, T) of the m-replica system for the molecu- 
lar liquid phase (note the crucial role of ergodicity in the 
liquid allowing for an unrestricted averaging over phase 
space; we cannot hope to do such a calculation for the 
amorphous solid with its restricted phase space). Pro- 
vided our system indeed follows the heuristics of an en- 
tropy crisis, we can cast this free energy density into the 
form Q with f*(m,T) the solution of the saddle-point 
equation JSJ . The critical replica parameter to* (T) then 
is determined by the condition /„ (m* (T) , T) — f m i n {T) 
or, cquivalently, ^[/^(m* (T)), T] = 0. Fortunately, we do 
not need to know a priori the form of the configurational 
entropy density S(f,T): calculating the derivatives of 
the free energy density F{m, T) with respect to to at 
to = to*(T), we easily find that 




{(T/m 2 )S[Mm,T),T] 
+ (dnfjdf. [/* - TS(/„, T)/m] } m= 



d m F(m,T)\ 



= 



FIG. 3: (a) Sketch of the free energy density F(m, T) for a in- 
fold replicated system as a function of the replica parameter 
m for different values of the temperature T. The smooth de- 
crease of the maximum in F(m) with increasing temperature 
T reflects the (weak) temperature dependence of / m i n (T). (b) 
Free energy density as a function of T after eliminating the 
replica parameter m. 



F(m, T) located at a value m,(T) > 1 (i.e., T c < T) tells 
us that the original system resides in the liquid phase. 
The true free energy density of the original system with 
m = 1 then is given by the value i<ii qu id = F{m = 1,T) 
assumed at the border of the interval. Finally, the tran- 
sition from one regime to the other at m. t {T c ) = 1 deter- 
mines the transition temperature T c . 



and 



d 2 m F(m,T)\ 



[(d m f*)/m] 



<0, 



C. Implementation of replica calculation 



where we have made use of the specific shape of the func- 
tion S(f,T), cf. Fig. El We thus conclude that the func- 
tion F(m,T) exhibits a maximum at to = to*, cf. Fig. 
13 and we can determine the critical replica parameter 
m*(T) directly from the free energy density F(m,T) in 
the molecular liquid phase without explicit knowledge of 
the entropy density S(f,T). 

In summary, we can adopt the following strategy in or- 
der to obtain a proper thermodynamic description of the 
glass transition associated with an entropy crisis: First, 
we have to compute the free energy density F(m,T) 
of the TO-fold replicated system in the molecular liquid 
phase and for arbitrary values of the parameter m. Next, 
this function has to be analytically continued to arbitrary 
(non-integer) values of to, in particular, values within the 
interval < to < 1. In a last step, we have to find the 
maximum of the function F(m, T) within the interval 
[0, 1]. If the maximum of the replica free energy F{m, T) 
is realized at m*(T) < 1 (i.e., T < T c , cf. Fig. 0), we 
conclude that the original system (with rn = 1) is in 
the glass phase and its free energy density is given by 



We consider a system of N identical particles, confined 
within a macroscopic box of volume V and described by 
the Hamiltonian if[xi,X2, ...,xat] = H[x]. The partition 
function for to uncoupled copies of this system reads 



(N\) m 1111 

' i=l a=l 



N m „ ,jj a m 

nn /^exp[-/?$> M 

-— 1 — 1 J a=l 



(9) 



-Fgiass = F[to*(T), T]. On the other hand, 



where (3 — l/T denotes the inverse temperature and s 
is the particle size (or the lattice spacing) which we set 
to unity hereafter, s = 1. Also, we define the density 
p = N/V, which we keep constant in the thermody- 
namic limit, as well as the mean interparticle distance 
L = p^ 1 / . The free energy density is given by the ex- 
pression F (to, 0) = —(\nZ m )/mf3N. Following the above 
strategy, we assume that the spatial positions of the par- 
ticles in different replicas are correlated, i.e., the particles 
arrange in 'replica molecules'; technically this is imple- 
mented via a transformation to center of mass (x^) and 
relative coordinates (u"), x° = x; + u°, and assuming 
a maximum in that the displacements u° remain bounded, |u"| L. In 
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addition, the coordinates have to satisfy the constraint 



E 1 

a=l 



(10) 



Finally, the assumption that this replica system resides 
in the molecular liquid state implies that the molecular 
positions x^ take on arbitrary values within the entire 
system volume. Hence, the replicated partition function 
can be written in the form 



Zm, — 



, DN 



Nl 



N 



N m 



1 V III n ilk 

i=la=l J 6=1 



i=l 



J[ / d D x % exp -/3^ff[x + u Q ] . (11) 



Here, the additional factor m DN is due to the change of 
variables xf = Xj + uf with the constraint (|10l) . 

Our problem now has assumed a form which is simi- 
lar to the replica representation of usual glass problems 
with quenched disorder. The molecular coordinates Xj 
play the role of the disorder parameters, while the dis- 
placements u" represent the dynamical variables. Our 
task then is to average over the disorder parameters {x^ } 
in order to arrive at the partition function Z m expressed 
through a new effective replica Hamiltonian H m [u a ], 



N m 



Z m = fl fl / d D M<} 5(jr uf) exp {-0H m [u°]) , (12) 

„-_i — i J b=l 



where the replica variables u° usually become coupled. 
The final integration over the dynamical variables uf 
then will provide us with the free energy density F(m, /?). 
Note that in order to assure the consistency of the cal- 
culation, we have to compute the mean squared dis- 
placement amplitude (u 2 ) 1 / 2 and verify that it is indeed 
smaller than the average distance between the particles, 
(u 2 ) 1 / 2 <C L. The above program then is carried out for 
the specific interaction introduced in the next section. 



calculations within the continuum limit. More specif- 
ically, we will treat the particles as point-like objects 
when integrating over positions. On the other hand, a 
proper hard core repulsion is required in order to inhibit 
the trivial collapse of the system into a pair of clusters 
separated by R with N/2 particles each. 

For later convenience we introduce the volume integral 



d D xU(\x\) 



R^<\x\<R 

V 



-R) 2 



R^ro 



(14) 



with Vq = 2roSo the volume of a spherical shell 
with radius R and width 2ro (here, Sd = 2-k d I 2 
R(D-V/T(D/2) ~ i?^" 1 ) is the area of the D-dimen- 
sional sphere with radius R). Finally, we will assume a 
particle density near to that of the liquid/crystal phase 
by fixing the mean particle separation L = p~ x l D close 
to the interaction radius R, L ~ R. 

Simple geometric considerations tell that in addition 
to the low-energy crystal configuration (characterized 
by space periodicity and a fixed number of nearest 
neighbours), the present model also develops numerous 
metastable low-energy states with a disordered/glassy ar- 
rangement of particles. Such states are inhomogeneous 
in space and exhibit a notably smaller average number of 
nearest neighbours as compared to the crystal. Also, geo- 
metric considerations tell that such glassy configurations 
are 'well separated' from the ordered crystal: their trans- 
formation into the crystal state would require a global (on 
the scale of the entire system) rearrangement of the par- 
ticles which would involve large energies. It is this type 
of random glass-like configurations which is at the focus 
of our further studies below. 



III. REPLICA FREE ENERGY 

We start from the expression (|llfl for the replica par- 
tition function in the form 



D. Model interaction 



We define the inter-particle interaction in the Hamilto- 
nian m in terms of the spherically symmetric potential 
U(|x|) (cf. Fig.HJ we define R ± =R ± r ) 



u(M)- 



R? 



o. 



i?_ < |x| < R + , 
|x| ?[R-,R + }. 



(13) 



We require the thickness 2tq of the spherical attractive 
shell to be small compared to the radius R, ro -c R. At 
the same time, we choose a shell thickness ro large com- 
pared to the particle size, which we choose equal to unity 
for simplicity, ro > 1; this allows us to perform all the 



Z,, 



DNirN N m 

m V — — 



Nl 



1\ lit n lib 



i—l a—1 



n 



V 



■ exp 



r, m IV 

{-EI 



6=1 



U[\(Xi-Xj 



a—1 



)+<-!]}, 



where = (u" — u°). For large N, we can approxi- 
mate Nl ~ N N and obtain the prefactor m DN V N /Nl ~ 
exp[DMn(m) - Mn(p)] with the density p = N/V. The 
replica partition function then takes the form 

N m r, rn 
i=\a=l J 6=1 

x(exp(-^[x;u])> x (16) 
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with the average over positions 

N 



($[x 1 ,...,x JV ]) x = JJ / -^$[xi,...,Xjv]. 



(17) 



The average of the exponential in Eq. can be calcu- 
lated with the help of a cumulant expansion, 



<exp(-/3i?[x;u])) j 



(18) 



exp 



-0(H), 



2! 



P 2 ((H 2 }}, 



3! 



«# 3 »> 



where ((H k )) x denotes the k-th cumulant of the Hamil- 
tonian; in particular, 

((tf 2 )) x ee (tf 2 } x - (H)l, (19) 
((tf 3 )) x ee {H 3 ) x -S(H 2 ) x (H) x + 2(H)l. (20) 

A. First- and second-order contributions 

Using the potential l|13|l with its volume integral (|14J) . 
we find the first-order term 

m jv 

= 2 E E ( u [K x * - x ;) + u £0>* (2i) 

a— 1 i,j=l 

= ~ m i — U ^^ ~ m - 

For the second-order cumulant (see Fig. 0[a) for a dia- 
grammatic representation), we find 

1 m AT 

«# 2 »* = i E E^ocxi-x^+u^o 



a,b=l i,j,k,l 

xf7[|(xfc-xj) + u: 



fci |J//Xi,Xi,Xj,x k - 



(22) 



Only correlated positions = (k,l) (i.e., connected 

graphs) give a finite contribution (since (U(\5xij + 
<.|)C/(|5x H + u^|)> = (f/(|x|))<[/(|x|)», hence 



* m N » 

<^ 2 »* = i E E/ 



d^x 

u 



a,fc=l i,j 

xC/(|x + u?.|) C^(|x + <|). (23) 

We assume that typical deviations { | u| 2 ) 1//2 are small 
compared to ro (this assumption has to be verified at 
the end of the calculation), allowing us to expand the 
potential U[\(x + u?-|] in the displacement u" ■, 

m Dm 

£ C/(l(x + u?-|) c mU(\x\) + £ U' a (M) E ««(«) 



D m 

-5 E ^(M) 



ij(/3)- 



(24) 



(a) 



(b) 



where o- 



-o = l/(lx+u°l) 



o X o 



where o X o =jlQ}x\) 



FIG. 4: Diagrammatic representation of the second-order cu- 
mulant; (a) only connected graphs give a finite contribution, 
(b) contributions after expansion in the displacement u" 3 -. 



Here, the indices a,/3 denote spatial vector components 
and U' a , U'^p are corresponding derivatives. The expan- 
sion (|24|l involves the small parameter u 2 /r 2 ,, allowing its 
termination at the second order, cf. Fig. Hlb) f° r a dia- 
grammatic representation. Its second term (linear in u) 
vanishes due to the constraint (|10fl and we arrive at the 
expression 



N 2 

((H 2 )) x ~— m 2 (C/ 2 (|x|)), 



(25) 



m N 



af3=l 



Using the definition of the potential t/(|x|), Eq. (|13|) . we 
obtain 



uUR-<M<R^ 



la/3 



1- 



R 



R 



(26) 



and else, where n a = x Q /|x| is the unit vector in the 
direction of x. Accounting for the smallness of ro, ro <c 
R, we find the positional averages 



M|X|) 2 )X =- /^[C/dxl)] 2 ^^, 



(27) 



<[/(|x|)[/^(|x|)) x = 



V 

V I n 



t/(|x|)C/^(|x|) (28) 
V 2 



V r 2 D 



and substituting these results back into Eq. (|25fl we ob- 
tain 



((H 2 )), 



NpV 2 

771 

2 



(29) 



EEE1 

i^j ct—1 a—1 



a ,3=1 



Going back to the original displacement coordinates u° 
and accounting for the restrictions (the second condition 
corresponds to a global shift of all particles in the system) 



7 



EciLi u t — an d u f = 0' we obtain the following 
contribution from the first- and second-order cumulants, 



-0(H), 



2! 



((H 2 )) x 



(30) 



NpV 
' 2 

' Drl 



D m N 



a — 1 a— 1 



Next, we concentrate on the three-point (loop) contribu- 
tion and first simplify the expression by redefining the 
coordinates , Xj , X& , 



TV 



(34) 



a= E E^+^-i] 

a,6,c— 1 i,j,k 

x£/[|x 2 + u^|]f/[| Xl - x 2 + u? fe |]) xlX2 . 

We expand in the displacement u, cf. Fig.[5Jc), and arrive 
at the expression 



B. Third-order cumulant 

Next, we determine the contributions from the third- 
order cumulant ((H 3 )) x = (S/2 + 2A) which contributes 
with two terms describing two-point (S) and three-point 
(A) correlations. These correspond to the connected dia- 
grams shown in Fig. Ef a) and can be written in the form 

m N 

5= E ^{U[\x+v§\]U[\x+i%\]U[\x+v%\]) x (31) 

a,b,c— 1 i,j 

and 

m N 

A= E E^[l^-^) + <,l] ( 32 ) 

a,6,c— 1 i,j,k 

xU[\(xj -Xfc) + u) fe !]f7[|(x fc -x,) + uy) XiXjXfc . 
Expanding the two-point contribution Eq. (|31|) in the 



3iV 



iV J m J (C/(|x 1 |)C/(|x 2 |)C/(|x 1 -x 2 |)) XlXa (35) 
E(^(|xi|)C/(|x 2 |)C/^(| Xl -x 2 |)) XlX2 



Q/3=l 



m AT 
<j=l ij 



Using the definition (|13|) of the potential and its deriva- 
tives (HUl, we find for the above disorder averages (again 
assuming that r$ <C R) 



([/(|x 1 |)[/(|x 2 |)C/(|x 1 -x 2 |)) Xl)X2 
([7(|x 1 |)C/(|x 2 |)C/^(|x 1 -x 2 |)) xli 



VqVi 
V 2 



(36) 
2 KM 



r 2 D V 2 



with Vi ~ R D ~ 2 rQ the intersection volume of two spher- 
ical shells with radius R and width r ; the three-point 
contribution then takes the form 



(a) 



(c) 



(b) o 



^ + 2 




-Q O- 

-3 — Q- 



K3 



O- 




-O O X O 

-o + 3 o 9 




+ 3 




A = -Np 2 V Vim 3 



(37) 



— n — m 



N 



EEE[<< 



a— 1 a— 1 i 



Replacing the density p with the distance parameter L 
and assuming that R < L, we find that the three-point 
contribution is reduced with respect to the two-point 
term (I33f) by the small factor 



pVt 



D R D ~ 2 r 2 



< 1. 



(38) 



k-t/^order cumulants 



FIG. 5: Diagrams contributing to third order in the interac- 
tion potential: (a) two-point and three-point (loop) diagrams, 
(b) contribution of the two-point diagram after expansion in 
the displacement u, (c) contribution of the three-point dia- 
gram after expansion in the displacement u. 

displacement u and integrating over x, we obtain up to 
second order in the displacement u (cf. Fig.^b)) 



Comparing the magnitude of the above diagrams, we 
note that all two-point terms (cf. Fig. EJ) appear with 
similar weight and we have to resum them; the two-point 
part of the k-th-ordei cumulant contributes with a term 
{-(3) k {{H k )) 2pt /k\ where (k > 2, see for the k = 1 
contribution) 



(-l) k ((H k )) 



2pt 



NpVo k 
m 



(39) 



-NpV m 3 



6 -E^m 2 
Dr 2 



D m N 



EEEK( Q) ]- (33) 



ot — l a—1 i 



Dr 2 



D m N 



'"TEEM 
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FIG. 6: Diagrammatic representation of the two-point con- 
tribution to the fc-f/i-order cumulant after expansion in the 
displacement u. 



which does not depend on the replica number m. When 
proceeding with l|43l) we have to keep in mind two re- 
strictions which are limiting its validity: i) our contin- 
uum approximation prevents us from reaching very low 
temperatures, and ii) the smallness of the mean square 
displacements (u 2 ) <C Tq requires a sufficiently large pa- 
rameter /3m or a sufficiently small temperature T. We 
proceed with a detailed analysis of these restrictions. 

D. Restrictions 



On the other hand, the loop diagrams remain small: 
a simple geometrical analysis shows that each k-th-order 
loop involves an additional small factor of order (k > 3) 



(p^o) fc " 3 (p^i) 



R\ D (ro 
L ) \R 

D(k-2) 



L 

R 



D 



ro 



< 1. (40) 



Thus, we conclude that for the present model with prop- 
erly chosen parameters ro <C R and R < L all loop-type 
contributions produce only small corrections as compared 
to the main terms arising from the two-point diagrams. 

It remains to sum up the two-point contributions, 
which corresponds to the substitution \(3m + (f3m) 2 /2! + 
(/3m) 3 /3! . . . — ► [exp(/3m) — I], and we obtain the parti- 
tion function averaged over disorder in the form 

1]} 
(41) 



Z m = cxv{DN\n{m)-N\Ti{p) + {NpVo/2)[e pm - 

N D s m p m 

>< nnln/^)<E«j 

6=1 

m 



i—l a—1 a—1 



x exp 



/3m 



Drl 



a=l 



A simple Gaussian integration over the displacements 



then results in the expression 



exp^ — — ln(m) 



JVh(p) 



N P y o r,flm 



ND(m-l) ln ipV Q p^ 0m 



Drl 



[e? m - 1] 



1] 



(42) 



and taking the logarithm we obtain (after rearrang- 
ing terms) the replica free energy density F(m, (3) = 
~{l/N(3m)\n{Z m ) in the form, 



F(m,f3) = F (P) 



pVo 

2/3ra 



[e? m - 1] 



(43) 



D 



In 



pV L 2 (3m 



Drl 



2/3m 

where we have separated the term 



(m — l)ln[e 



[3m 



^ }' 



FoW - gin 



pVo(3 
Dr 2 



(44) 



In order to justify the continuum limit in the integra- 
tions over the displacements u% a ^ , cf . Eq. 141J1 , where we 
assume that |n°( a )| ^ 1, we have to demand that the 
width 



Drl [( 



/3 m 



-1]<C1. 



(45) 



Hence, the application of the continuuum limit puts a 
lower limit on the system temperature, or, alternatively, 
on the largeness of the replica parameter m. We will 
return to this point in the next section. 

The second restriction derives from the condition that 
typical values of the displacements as defined by 

the Boltzmann weight in Eq. (|41|l should be small on the 
scale of the width ro of the attractive shell; the Gaussian 
integration then yields the condition 



m * — ' 

a—1 



2 ) 



Dr 2 /pVo 



/3m[exp(/3m) — 1] 



<-rl. 



Since pV ~ {R/L) u (r /R) < 1 for R < L, we can 
conclude that the result for the replica free energy density 
describing the molecular liquid state is valid if (3m is 
sufficiently large and bounded by the condition 



pVo /3m exp (/3m) ^> 1. 



(46) 



The violation of this condition implies that the particles, 
originally assumed to be bounded in 'replica molecules', 
would escape from the potential well ?7(|x|) and become 
effectively decoupled. In this situation the state of our 
replica system would not correspond to the molecular 
liquid phase any more and the above analysis cannot 
be applied. The same conclusion also follows from our 
general qualitative arguments in Sec. II. B: the limit of 
small replica parameter m corresponds to a high effective 
temperature of the replicated system where no 'replica 
molecules' could survive. 



IV. RESULTS AND CONCLUSIONS 
A. Free energy density and T c 

Let us return to the free energy density l|43|l and 
determine the replica parameter m. Assuming that 
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exp ((3m) ^> 1, we first simplify the function F(m,(3), 



F(m,/3) ~ F (/3) 



e /3m 



-In 



pV L 2 (3m 



2(3m 2/3m 
+D(m -l)/2, 



(47) 



and find its maximum in the replica parameter m from 
the condition d m F(m, 0) — 0, 



^(/3m-l)e' 3m 
D y ' 



(pr. 



- 1 = In 



pV a L 2 (3m 



Drl 



(48) 



With the parameters r /L, pVo <C 1 and pVoL 2 /Vg ^> 1, 
the above equation assumes the solution (to logarithmic 
accuracy) 



0m J (3) ~ In— +ln 
pV 



\n{pV L 2 /eDr 



HD/pVo) 



1 ^ D 



(49) 



In the following, we concentrate on the case of matching 
density p ~ R~ D , i.e., L ^ R, and substituting p ~ R~ D , 
Vq ~ R D ^ 1 ro we easily verify that the above conditions 
are satisfied, 



pV 
L 2 



D 



^ ~ 1° « 1 
F J? 



'•() 



F 
'•() 



> 1. 



The solution m* then assumes the simplified form 
DR 



0m*(l3) -ln- 



ln 1 



ro 



/? r J 



(50) 
(51) 

(52) 



The transition temperature T c is defined by the condition 
m*(/3 c ) = 1 (see Sec. II. B) and we arrive at the estimate, 



1 



HDR/r Q ] 



< 1. 



(53) 



Next, we should verify that the above conditions l|45|l and 
(1461) are satisfied; they can be cast into the form 



with 



1 < — /W 3 " 1 < mDrl; 
R 



(3m > /3m* (T) - ln(L>F/r ). 



(54) 



(55) 



The first relation (guaranteeing a small displacement am- 
plitude (u 2 ) <C r-p) is satisfied within the entire low- 
temperature phase < T < T c . The second relation 
(allowing us to make use of the continuum limit) implies 
that T 1/tq. This condition is satisfied at T c since 
ln(DR/ro) <C Tq for sufficiently large ro; however, it is 
violated at low temperatures T — 1/Vq , thus limiting the 
applicability of our results to sufficiently high values of 
T. 

The logarithmic dependence of T c found above can be 
easily understood from an order of magnitude estimate 



of the free energies of the liquid and glass phase. Let us 
assume that the average number of particles interacting 
in the frozen state is of the order of D (see below, Eqs. 
(|73[l and 174|) ). Then each particle contributes with an 
energy —D to the free energy of the system. Second, 
the freezing into the glass state confines the particle to 
the volume which is small with respect to the volume 
R D available in the liquid state. This produces a deficit 
\n(R D /tq ) in the entropy of the frozen state as compared 
to the liquid. The difference in the free energies between 
the frozen and liquid states then can be estimated as 
AF — — D + T\n(R D fr®); this quantity turns negative 
at temperatures T < T c — [ln(F/ro)] _1 and the frozen 
state becomes preferable. 

An independent confirmation of the result Ij53(l can be 
obtained from a virial expansion^: to third order the 
equation of state takes the form 



pV 
NT 



1 + B(T)^- + C{T) 



V 



N 2 

] T 2 



(56) 



with the coefficients 



B{T) = ~\J d 3 xf(x), (57) 
C(T) = -\j rf 3 xid 3 x 2 /(x 1 )/(x 2 )/(x 1 -x 2 ), 

and /(x) = exp[— (3U(x)] — 1. Inserting the expression 
(|13[1 for the potential U we can rewrite 



/(x) 



exp(/3), R_ < |x| < F + , 
0, |x|£ [F_,F + ], 



(58) 



and using the results (fT4*jl and (|3l)|l above, we find B ~ 
-R D {r /R)exp{(3) and C(T) ~ -R 2D {r /R) 3 exp(3/3). 
Assuming a high particle density with L ~ F, we find 
that the virial expansion diverges when (tq/R) exp(/3) — 
1, thus reproducing the critical temperature lj53|) . 



B. Configurational entropy 

We can use our free energy expression F(m, (3) to con- 
struct the form of the configurational entropy S(f,T), 
cf. Fig. 1. Taking the derivative of the free en- 
ergy F(m, T) = f(m,T) - (T/m)S[f{m,T),T] and of 
mF(m, T) with respect to m and using that dfS = m/T, 
cf. and 0, we obtain the relations 



—d m F(m,T) = S(m,T), 
a m [mF(m,T)]=f(m,T), 



(59) 
(60) 



and eliminating the parameter m, we arrive at an ex- 
pression for S(f,T). We use the free energy lj4l?)l in the 
limit of large (3m and expand (to third order in m — m*) 
around the maximum F(m*(T), T). Expressing m — m» 
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through /- / min (with / min = F(m*(T),T)) we find that 
the configurational entropy near / m i n assumes the form 

S(f, T) a /3m* (/ - / min ) - ^L"^* (/ - / min ) 2 . (61) 

This result exhibits the shape expected for a configura- 
tional entropy triggering an entropy crisis as discussed in 
Sec. II.B, see Fig.0 



C. Order parameter 

The above formal replica computations have been per- 
formed along the lines outlined in section II and fit well 
the original physical ideas regarding the nature of the 
liquid-glass phase transition. In order to describe the 
phase transition in more quantitative terms one has to in- 
troduce a properly defined order parameter which should 
be a measureable quantity e.g., in computer simulations. 
Such an order parameter is easily defined if we exploit 
the replication trick where we have introduced m iden- 
tical copies of the same system. Let us assume that all 
these systems are allowed to thermalize independently 
but with the same (random) starting positions of the 
particles. If, at sufficiently high temperatures, the ther- 
modynamic state of the system is a liquid, we expect the 
spatial positions of the particles belonging to different 
systems to be uncorrelated. On the other hand, if the 
thermodynamic state of the system is a frozen glass with 
each particle localized in a limited part of space, then the 
positions of the particles belonging to different replicas 
remain correlated, although the systems are uncoupled. 
Keeping in mind this qualitative scenario, we introduce 
the correlator 



Changing variables according to x° = Xj + u", we find 



1 N 
G = ~Y( 



<.2) 



(62) 



where x^ 1 ' is the position of the i-th particle in the system 
number 1, and x^ 2) is the corresponding position (again 
of the i-th particle) of the system number 2. If the ther- 
modynamic state of the system is a liquid, then the value 

of G is proportional to TV -1 J2iLi(( x 'i 1} ) )j which is of 
the order of the linear size squared of the system and 
becomes infinite in the thermodynamic limit. 

On the other hand, if the system is in the glass state, 
the situation is very different. In this case both x^ 1 ' and 
Xj are localised near the same equilibrium positions and 
G remains finite. Since all m systems (replicas) consid- 
ered here are equivalent, it is convenient to write the 
correlator l|62|l in the symmetric way 



G 



1 N 

-Y 



— J m(m — 1) 



) 2 ). 



.,6=1 



We compute G in the glass phase following the procedure 
described in section II and implemented in section III. 



G = 



m(m— 1) 



a, 6=1 

2Dm „. a . 



(63) 



with the distribution of displacements it" given by Eq. 
(|41|l . A simple Gaussian integration yields 



pVo 



(64) 



where pVo ~ r /R (see Eq. 1(50)1) and the value of 
m^T) is given by the saddle-point solution /3m* (T) ~ 
In (DR/ro). Thus, in the glassy phase we find the finite 
value 



G ~ Tr 2 . 



(65) 



We then define the order parameter Q — (1 + G) _1 as- 
suming a finite value in the glassy phase and vanishing 
in the liquid, 



0. 



in the liquid, at T > T c 



Q = 



1 (66) 
2, in the glass, at T < T c , 



l + Tr z Q 



where the transition temperature T c is given by Eq. H53|) . 
Note that the value of Q remains finite at T — T c , as 
expected for a phase transition driven by an entropy cri- 
sis. On the one hand, the physical order parameter de- 
scribing this phase transition exhibits a finite jump at T c 
(as expected for a first-order phase transition), while, on 
the other hand, the free energy of the system is continu- 
ous in the transition point (with a continuous derivative 
as in a second-order transition). On a qualitative level 
this is easily understood: in the present approach the 
glass phase is characterized by the typical size of spatial 
cells where the particles remain localized. This volume 
is small in the low-temperature limit and grows with in- 
creasing temperature. At the transition to the liquid, the 
localization length is of the order of the interparticle dis- 
tance and thus remains finite. Beyond the transition the 
particles move freely in the liquid phase and the local- 
ization length is infinite. Hence the transition resembles 
the usual solid-liquid transition, however, with the jump 
in entropy (latent heat) replaced by the entropy crisis 
scenario, guaranteeing the smooth transition in the free 
energy density. 



D. Free energy and entropy densities 

Finally, we analyze in some more detail the value of the 
free energy- and entropy densities of the liquid and glassy 
phases. The free energy density F{m, f3) of the liquid 
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derives from the result l|43(l with the replica parameter 
m = 1 (we assume that (J ^> 1 and choose L ~ R), 



-F|iquid(/3) 



TO 

2R(3' 



the entropy density is given by the derivative 

Sliquid iP) — /^ 2 d/3-fliquid 



(67) 



(68) 



We observe that at sufficiently low temperatures the en- 
tropy density of the liquid would become negative. For- 
mally, this takes place at j3 > f3 c with f3 c defined by the 
condition 



(/3 C 



ro 



(69) 



This negative entropy then signals the presence of the 
glass (or entropy crisis) transition: with Vq/R <C 1 we 
obtain an estimate for the transition temperature, 



the prefactor A assumes a value close to unity, 

InHDR/rp) 
~ 2ln(DR/r ) ' 



(74) 



We thus find that our glass phase is 'well packed' with 
slightly more than D particles in a (interacting) near- 
est neighbor position on average. Note, that the entropy 
l|72p becomes negative at sufficiently low temperatures, 
i.e., when ~ r$ ^> 1. This unphysical result is a con- 
sequence of the breakdown of our continuum approxima- 
tion at these low temperatures, cf. Eq. (|45|l . 

Next, we analyze the situation near the glass transition 
temperature. A general expression for the entropy of the 
glassy phase (without assuming that T <C T c ) can be 
obtained from l|47|l . 



Sglass(/3) 



g, p 

2 n DRr Q 
2R 



D\nR 



(75) 



D 



lnm»(/3), 



R 1 DR 

f3 c ~ In , 

ro 



(70) 



which is in agreement with the analytic result (|5I 

The above arguments do not imply that the entropy 
of the liquid at the phase transition is equal to zero. In 
fact, as we shall see below, the phase transition into the 
glassy phase takes place before the entropy becomes zero. 
The above consideration of the liquid entropy is merely a 
qualitative estimate of the temperature below which we 
would run into trouble if we were to use the result (|67|) . 

Turning to the frozen state, we first discuss the situ- 
ation deep in the glassy phase. The free energy density 
at (3 ^> (3 C follows from the result l|47|) with m — m*((3) 
defined by (J52J|. Assuming L ~ R we find the simplified 
expression 



F.. 



glass 



D 
~2 



i _ In Mn{DR/r ) 



0) ]n(DR/r ) 



1 r 2 
--ln^° 

P P. 



(71) 



Taking the derivative with respect to temperature we ob- 
tain the entropy density of the glassy phase, 



S g la SS (/3) = p 2 d P F gl . 

D 
~ ~2 



(72) 



1 



ln^- 

P 



Mn(DR/r ) 



HDR/ro) 



Combining the results for the free energy and the en- 
tropy, we can derive an expression for the average energy 
per particle in the glassy phase which scales with the 
dimension D, 



E = F< 



glass 



glass 



AD: 



(73) 



where we have used the relation (|48|l in order to allow for 
a simple comparison with the result l|68l) of the liquid. At 
T c we have m* = 1 and the entropy of the glass coincides 
with that of the liquid, 



Sglass(Pc) = ^liquid (Pc) 



(76) 



= D\nR-^{p c -l)e^. 



Inserting the estimate for f3 c , Eq. (|53|l . we find that the 
entropy is positive at the transition point, 



S s iass{Pc) = <S'ii qu id(/3c) ~ -Dlni? - 



D i R n 
—In— > 0. 

2 r 



(77) 



Another remark concerns the low-density or gas limit 
of our model away from the densely packed limit with 
L ~ R considered above. The system then does not 
exhibit any signature of the above phase transition. For- 
mally, we find no maximum in the replica free energy 
density (|43|l : rather, F(m,[3) diverges for m — * 0. How- 
ever, in this limit the constraint l|46() is violated and our 
calculation is not valid any more. In fact, the constraint 
(|46|l guarantees the smallness of the thermal displace- 
ment amplitude (u 2 ) < Tq; its violation tells us that our 
basic assumption of a molecular liquid phase is corrupted, 
which signals that the original system is in the liquid 
state. 



E. Conclusion 

In conclusion, we have derived analytic results for the 
liquid-glass phase transition in a model glass former, us- 
ing the replica mean-field theory proposed by Mezard and 
Paris: 1 -! based on the entropy crisis scenario. Of course, 
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the present study does not pretend to describe the physics 
in actual realistic glasses, which probably is much more 
complicated and furthermore essentially non-equilibrium 
in nature. Nevertheless, our analysis demonstrates that 
there exists a class of statistical systems which, on the one 
hand, incorporate essential features of structural glasses 
and, on the other hand, can be studied within an equilib- 
rium statistical mechanics approach, similar to numerous 
other disordered systems. Although the present calcula- 
tion was limited to a mean-field approximation, we have 
obtained a physically consistent set of results: Our repli- 
cated system indeed exhibits a free energy F(m, T) with 
a maximum m*(T) within the interval [0, 1] at low tem- 
peratures. This maximum shifts towards unity with in- 
creasing T and the condition m*(T c ) = 1 defines a glass 
temperature T c ~ l/\n(DR/ro) which is in agreement 
with various estimates. The maximum disappears upon 
diluting the system, indicating the persistence of the liq- 
uid state for low densities p <C R~ D . The behavior of the 
configurational entropy S(f, T) is fully consistent with 
the freezing scenario based on an entropy crisis and we 
recover the expected characteristics of a smooth transi- 
tion in the free energy combined with a jump in the order 
parameter. In addition, our analysis provides a struc- 



tural information on the glass state: calculating the mean 
particle energy from the free energy and entropy expres- 
sions, we find that the system freezes into a glass state 
with particles binding to slightly more than D neighbors. 
All these results have been derived in analytic form, pro- 
viding direct access to the parametric dependence of the 
physical results. 

On the other hand, we have to admit that at the 
present stage it is difficult to juge which of the features 
of the phase transition scenario found here are specific 
to the particular model and its mean-field solution and 
which of them are generic and reflect the general situation 
encountered in the structural glass transition. At least 
part of this proviso can be handled by pushing the anal- 
ysis beyond the mean-field approximation. With the ex- 
perience gained in the present study, we believe that the 
theoretical construction of this next step, which would 
include the effect of thermodynamic fluctuations, does 
not look unrealistic. 

We thank Lev Ioffe and Serguei Korshunov for interest- 
ing discussions and acknowledge financial support from 
the Swiss National Foundation through the Center of 
Theoretical Studies at ETH Zurich. 
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